Low-temperature magnetization of (Ga,Mn)As semiconductors 
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We report on a comprehensive study of the ferromagnetic moment per Mn atom in (Ga,Mn)As 
ferromagnetic semiconductors. Theoretical discussion is based on microscopic calculations and on 
an effective model of Mn local moments antiferromagnetically coupled to valence band hole spins. 
The validity of the effective model over the range of doping studied is assessed by comparing with 
microscopic tight-binding/coherent-potential approximation calculations. Using the virtual crys- 
tal k • p model for hole states, we evaluate the zero-temperature mean-field contributions to the 
magnetization from the hole kinetic and exchange energies, and magnetization suppression due 
to quantum fluctuations of Mn moment orientations around their mean-field ground state values. 
Experimental low-temperature ferromagnetic moments per Mn are obtained by superconducting 
quantum interference device and x-ray magnetic circular dichroism measurements in a series of 
(Ga,Mn)As semiconductors with nominal Mn doping ranging from ~2% to 8%. Hall measurements 
in as-grown and annealed samples are used to estimate the number of uncompensated substitutional 
Mn moments. Based on our comparison between experiment and theory we conclude that all these 
Mn moments in high quality (Ga,Mn)As materials have nearly parallel ground state alignment. 



PACS numbers: 75.50.Pp,75.30.Gw,73.61.Ey 

I. INTRODUCTION 

Early experimental studies of (Ga,Mn)As ferromag- 
netic semiconductors, reporting large apparent magneti- 
zation deficits?^ motivated a theoretical search for pos- 
sible intrinsic origins of frustrating magnetic interactions 
in this material. Using a wide spectrum of computational 
techniques, ranging from ab initio LDA methods 3 ^ and 
microscopic tight-binding approximations 5 " to semiphe- 
nomenological, k ■ p kinetic-exchange models^i2i2ii£ the 
theoretical studies have identified several mechanisms 
that can lead to non-collinear ground states. The ob- 
servation that long wavelength spin-waves with negative 
energies frequently occur within a spherically symmetric 
kinetic-exchange model illustrates^ that randomness in 
the distribution of Mn moments can result in an instabil- 
ity of the collinear ferromagnetic state. Frustration can 
be further enhanced when positional disorder is combined 
with anisotropies in Mn-Mn interactions. The pd char- 
acter of electronic states forming the magnetic moment 
leads to magnetic interaction anisotropies with respect 
to the crystallographic orientation of the vector connect- 
ing two Mn moments When spin-orbit coupling is 
taken into account magnetic interactions also be- 
come anisotropic with respect to the relative orientation 
of the Mn-Mn connecting vector and the magnetic mo- 
ment. 

Some degree of non-collinearity is inevitable as a com- 



bined consequence of positional disorder and spin-orbit 
coupling. Nevertheless a large suppression of the ferro- 
magnetic moment is not expected theoretically 5 in metal- 
lic (Ga,Mn)As samples with Mn concentrations above 
1%. The minor role of non-collinearity is due largely to 
the long-range character of magnetic interactions, which 
tends to average out the frustrating effect of anisotropic 
coupling between randomly distributed Mn impurities. 
In this paper we present detailed calculations of zero 
temperature magnetization in (Ga,Mn)As ferromagnets 
and compare the results with superconducting quantum 
interference device (SQUID) and x-ray magnetic circu- 
lar dichroism (XMCD) measurements in a series of sam- 
ples with nominal Mn doping ranging from ~2% to 8%. 
Our calculations neglect effects that would lead to non- 
collincarity, appealing to expectations that these effects 
are small. Our assumption is consistent with experi- 
mental observations of larger ferromagnetic moments in 
recently synthesized high-quality samplesii*^ The sub- 
stantial magnetization suppression seen in many early 
(Ga,Mn)As samples is attributed here primarily to the 
role played in those samples by interstitial Mn atoms. 
The consistency of the theoretical and experimental data 
that we are able to achieve, allows us to rule out any 
marked magnetic frustrations in the ground state of high- 
quality (Ga,Mn)As ferromagnets, and helps to clarify the 
character of magnetic interactions in this material. 

Two distinct theoretical approaches are used in the pa- 
per to discuss magnetization in (Ga,Mn)As semiconduc- 
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tors. In the more microscopic approach we account ex- 
plicitly for the five d-orbital electrons on a substitutional 
MnQa impurity, and for the strong on-site Coulomb cor- 
relations that suppress spin and charge fluctuations of 
the L — 0, S — 5/2 state of the atomic Mn d-shell. 
Magnetism in the mixed crystal arises in this picture 
from electron hopping between the Mn d-states and p- 
orbitals concentrated on the As sublattice that form the 
top of the host semiconductor valence band. For weak 
p — d hybridization, a second approach is possible. The 
Schrieffcr- Wolff transformation 13 allows us to map the 
microscopic Hamiltonian onto an effective Hamiltonian 
for local S = 5/2 moments and valence band states 
whose coupling is described by the kinetic-exchange term, 
JpdS(R) • s(r)<5(R — r), where S and s are the local mo- 
ment and valence band state spin operators, respectively. 
This approach will fail if the p — d hybridization is too 
strong, but appears to be reliable for (Ga,Mn)As. 

The paper is organized as follows: In Section [n] 
we identify the key physical considerations related to 
ground-state magnetization of (Ga,Mn)As ferromagnets 
by focusing first on a single Mn(d 5 +hole) complex. We 
recall the connection between p — d hybridization and 
the antiferromagnetic kinetic-exchange coupling, digress 
on the sign of the hole contribution to total moment per 
Mn, and discuss the expected mean-field contribution to 
magnetization per Mn from the Mn local moment and 
from the antiferromagnetically coupled hole. We also ex- 
plain that quantum fluctuations around the mean-field 
ground state are generically present because of antiferro- 
magnetic character of the p — d kinetic exchange interac- 
tion. 

Magnetization calculations for the many-Mn-impurity 
system are discussed in Section IIIII The relevant 
considerations here parallel those that apply for iso- 
lated Mn(d 5 +hole) complexes, but differ in detail be- 
cause of interactions between moments. Zero temper- 
ature ferromagnetic moments per Mn are first studied 
within the tight-binding/coherent-potential approxima- 
tion (TB A/CPA) model. The results of these microscopic 
calculations indicate that (Ga,Mn)As is in a weak p — d 
hybridization regime over the whole range of Mn concen- 
trations that we study. These calculations help establish 
theoretically the validity of the effective kinetic- exchange 
model. The virtual crystal approximation and the k • p 
effective Hamiltonian are then used to evaluate contribu- 
tions to the mean-field magnetization from hole kinetic 
and exchange energies and to confirm the expected weak 
role quantum fluctuations around the mean-field many- 
body ground state. 

Experimental SQUID and XMCD data are presented 
in Section IIVI Partial concentrations of substitutional 
and interstitial Mn impurities and the corresponding 
number of uncompensated local moments are derived 
from the nominal Mn doping and from Hall measure- 
ments of the hole density in as-grown and annealed 
samples^ The collinearity of the ferromagnetic ground 
state in the (Ga,Mn)As materials we study is tested by 



comparing experimental data with theoretical calcula- 
tions. Section briefly summarizes main conclusions 
of the paper. 



II. MAGNETIZATION OF AN ISOLATED 
Mn(d 5 +HOLE) COMPLEX 

A. d 5 +hole picture 

Most of the spectral weight near the top of the GaAs 
valence band originates from As 4p levels. Magnetic cou- 
pling between these states and the strongly localized 3d 
electrons on the substitutional MnQ a atom, which is at 
the origin of ferromagnetism in (Ga,Mn)As materials, is 
dominated by the p — d hybridization contribution^ The 
majority spin d-shell, with all five electron spins aligned 
forming a large spin S = 5/2, has its spectral weight cen- 
tered ~ 3.5 eV below the top of the valence band. The 
empty, minority-spin c?-level overlaps with the conduction 
band. Hybridization between p and d orbitals therefore, 
increases the energy of the majority-spin p orbitals (level 
repulsion) and favors the occupation of minority-spin p 
orbitals. This is the basic origin of the antiferromagnetic 
interaction represented by J pc i, as explained further in 

Fig.m 

Apart from providing a magnetic moment, the sub- 
stitutional MiiQa impurity acts as a moderately shallow 
single-acceptor in GaAs^ The cartoon in top panel of 
Fig. shows the splitting of the impurity level due to 
p — d hybridization and the lower panel of Fig. ^ illus- 
trates the many-Mn system in which the impurity level 
is broadened, eventually merging with the host valence 
band at higher Mn concentrations. To avoid confusion 
that may result from using the hole picture to describe 
magnetization of carriers in p-type (Ga,Mn)As materi- 
als, we make a digression here and explain the relation 
between magnetizations as evaluated using the physical 
electron-picture (as in Fig. ^) and magnetizations evalu- 
ated using the indirect but computationally more conve- 
nient hole-picture. 

Magnetization at T — is defined thermodynamically 
by the dependence of the ground-state energy E on ex- 
ternal magnetic field B: 



In this paper we always assume B || +z. In a mean- 
field picture the magnetization is related to the change 
of single-particle energy with field, summed over all oc- 
cupied orbitals. Orbitals that decrease in energy with 
field make a positive contribution to the magnetization. 
For B || + z, the d-electron spins are aligned along — z- 
direction (down-spins) and the majority spin band elec- 
trons have spin-up due to antiferromagnetic p — d ex- 
change coupling. Then, if the majority spin band moves 
up in energy with B and the minority band moves down, 
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as illustrated in the left part of Fig. [3 the band kinetic 
energy increases with B and, according to Eq. (JIJ, the 
corresponding contribution to the magnetization is neg- 
ative. In the hole-picture, we obtain the same respec- 
tive sense of shifts of the majority hole and minority hole 
bands, as illustrated in the right part of Fig. and there- 
fore the correct (negative in our case) sign of the magne- 
tization. The cartoon shows that in order to circumvent 
the rather confusing notion of a spin of a hole in magneti- 
zation calculations, the full Hamiltonian H{B) should be 
derived in the physical electron picture where the sign of 
the coupling of the electron spin to B is unambiguously 
defined. The electron picture — > hole picture transforma- 
tion (H(B) — > —H{B)) and the clearly defined notion of 
majority and minority bands in either picture guarantees 
the consistency in sign of the calculated magnetization. 
Note that the language used here neglects spin-orbit in- 
teractions which lead to single-particle orbitals that do 
not have definite spin character. Although spin-orbit in- 
teractions are important they can be neglected in most 
qualitative considerations, like the ones we explain here. 
In this paper we occasionally make statements which ne- 
glect spin-orbit interactions, and they should always be 
understood in this spirit. 

The electron-electron exchange energy has a nega- 
tive sign and its magnitude increases monotonically 
when moving from the paramagnetic to the half-metallic 
(empty minority band) state. This together with Eq. 
implies that the magnetization contribution from the 
electron-electron exchange energy has the same sign as 
the contribution from the kinetic energy. Using the same 
arguments as above we see in the exchange energy case 
the sign of magnetization is also treated consistently in 
the electron picture — > hole picture transformation. 
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FIG. 1: Electron-picture cartoon: splitting of the isolated Mn 
acceptor level (top panel) and of the top of the valence-band 
in the many-Mn system (bottom panel) due to p-d hybridiza- 
tion. 




FIG. 2: Cartoon of Zeeman coupling of an external magnetic 
field assuming g > in the electron and hole pictures for our 
valence band coupled to Mn moments system. Majority band 
in both electron and hole pictures moves up in energy result- 
ing in a negative band-contribution to the magnetization. 

B. Mean-field magnetization 

In a model which for the many-Mn-impurity system 
corresponds to a mean-field approximation, the ground 
state wavefunction of the Mn(d 5 +hole) complex reads 
\S Z = —S)\j x = +j) and the magnetization per Mn 
equals rriMF = (gsS — gjj)^B, where S and j are local 
d-electrons and hole moments and gs and gj are the re- 
spective Lande g-factors. The five d-electrons have zero 
total orbital angular momentum, i.e. gs — 2, and for 
spin j — 1/2 hole (gj = 2) we get tumf = 4/zb. Hole 
states near the valence band edge have a p-character, 
however, so more realistically we should consider gjj = 
4/3 * 3/2 = 2 which gives vumf = 3/jb. We show in the 
next section that this basic picture of a suppressed tximf 
due to holes applies also to highly Mn-doped (Ga,Mn)As 
materials although the magnitude of the mean-field hole 
contribution is weaker because of the occupation of both 
majority and minority hole bands and, partly, because of 
spin-orbit coupling effects. 

C. Quantum fluctuations contribution to 
magnetization 

The two-spin, S and j model allows us also to read- 
ily demonstrate the presence of quantum fluctuations 
around the mean-field ground state, which is related to 
the antiferromagnetic sign of the S • j coupling. We show 
that for the Mn(d 5 +hole) complex, quantum fluctuations 
are expected to weakly suppress the mean-field magneti- 
zation 77i mf- Detailed many-body calculations discussed 
in the following section confirm the role of quantum fluc- 
tuations is also weak for the many-Mn systems. 

In the limit of B — > we can write the two-spin Hamil- 
tonian as, 

H = JS.j=±(S? ot -S 2 -3 2 ), (2) 
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where Stot = S + j. For comparison we first consider 
ferromagnetic coupling, J < 0. Since in this case S to t = 
S + j, the ground-state eigenenergy, Efm — — ^[(S + 
s)(S + J + 1)- S(S + 1) - j(j + 1)] = -\J\Sj, equals to 
the mean- field energy, i.e., the mean- field state is exact. 
For ferromagnetic S • j coupling, the quantum fluctuation 
contribution to the magnetization is strictly absent in the 
many-Mn case only when the hole system is half-metallic 
(i.e. when the minority band is empty). We can see this 
by introducing spin raising and lowering operators in the 
Hamiltonian J5J, 



H = J 



s z ] z + -(s + r + s-f 



(3) 



Quantum fluctuations are absent when the transverse 
spin terms above annihilate the many-particle ground 
state. When acting on a state with all localized spins and 
all band spins polarized in the same direction, both trans- 
verse terms produce zero. For partially spin-polarized 
bands, quantum fluctuation corrections, although not 
strictly zero, are qualitatively smaller than in the an- 
tiferromagnetic case. 

For antiferromagnetic coupling (J > 0), Stot = S — 
j and the corresponding ground-state energy Eaf = 
l 4[(S-j)(S-j+l)-S(S+l)-j(j+l)} = -\J\(Sj+j) is 
lower than the mean-field energy. The mean-field ground 
state is not exact here and quantum fluctuation correc- 
tions to the magnetization will be non-zero in general. To 
estimate the correction we write the exact ground-state 
wavefunction as, 

\i>) = \s to t = s - j,s to t,z = -(s - j)) 



s 



S + J 



-\S Z = -S,j z = +j) 



-\S, = -S + l,j, = +j-l), (4) 



V S + J 

the mean-field wavefunction as, \i/))mf — \S Z — —S,j z — 
+j), and evaluate the respective expectation values of 
the Zeeman Hamiltonian, gs^sBS z + gj[iBBj z . From 
Eq. Q we then obtain that the difference between the 
exact and mean-field state magnetizations is given by 



m - m M F = m QF 



Ms 



J 



S + j 



(9s-9j) 



(5) 



When j — 1/2 and gs — gj — 2 the quantum fluctuation 
correction to the magnetization vanishes even though 
the mean-field ground state is not exact. The quantum 
fluctuation correction to the magnetization remains rel- 
atively weak also even when we adopt the more realistic 
description of the valence-band hole moment (j = 3/2, 
gj = 4/3), for which itiqf = — 0.25/is. 



III. MAGNETIZATION OF 
MANY-Mn-IMPURITY SYSTEM 



A. Moments per Mn in the microscopic TB A/CPA 

model 



The TBA description of (Ga,Mn)As mixed crystals is 
particularly useful for explaining the complementary role 
of local and itinerant moments in this p-type magnetic 
semiconductor. The language that is used to describe this 
interplay can differ depending on whether a fully micro- 
scopic or a kinetic-exchange model is employed, and this 
difference has sometimes led to confusion. This section 
represents an attempt at clarity. At the same time we 
find that the TBA/CPA results help establish the valid- 
ity of the antiferromagnetic p — d kinetic-exchange model 
and of the virtual crystal approximation for describ- 
ing collinear ground states in highly doped (Ga,Mn)As 
ferromagnets ,A2ii& 

In the TBA/CPA calculations, the hole density is var- 
ied independently of Mn doping by adding non-magnetic 
donors (Si or Se) or acceptors (C or Be). The d-electron 
magnetic moments of all Mn atoms are aligned along +z- 
axis. The parameterization of the TBA Hamiltonian was 
chosen to provide the correct band gap for a pure GaAs 
crystal^ and the appropriate exchange splitting of the 
Mn d-states. Local changes of the crystal potential at 
Mn, represented by shifted atomic levels, were estimated 
using Ref . 0. Spin-orbit coupling is not included in our 
TBA Hamiltonian. In the CPA, disorder effects appear 
in the finite spectral width of hole quasiparticle states. 
The TBA/CPA technique can, therefore, capture changes 
in the p — d interaction with doping due to both chemi- 
cal alloying effects and positional disorder. In Fig. [3] we 
show the microscopic TBA/CPA magnetic moments per 
Mn, ttitba, in Gai-^Mn^As ferromagnets plotted as a 
function of the hole density p relative to the Mn con- 
centration A^Mn = 4x/aj i c (ai c is the semiconductor host 
lattice constant). The ttitba is obtained here using the 
electron picture by integrating over occupied states up to 
the Fermi energy. 

A common way of microscopically separating contribu- 
tions from local atomic and itinerant moments is by pro- 
jecting the occupied electron states onto Mn d-orbitals 
and sp-orbitals, respectively. In this decomposition, the 
resulting local Mn moments are smaller than 5/is per Mn 
due to the admixture of d-character in empty states near 
the valence band edge. The effective kinetic-exchange 
model employed in the following sections corresponds, 
however, to a different decomposition of contributions, in 
effect associating one spectral region with local Mn mo- 
ments and a different spectral region with itinerant hole 
moments. The kinetic-exchange model, in which local 
moments have S — 5/2, is obtained from the microscopic 
TBA/CPA model by expressing the total TBA/CPA mo- 
ment as the difference between a contribution rn™ BA re- 
sulting from integrating over all electronic states up to 
mid-gap, i.e. including the entire valence band, and a 
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contribution corresponding to the integral from Fermi en- 
ergy to mid-gap. As long as the valence-conduction band 
gap is non-zero, the former contribution is independent 
of valence band filling and equals to the moment of an 
isolated Mn atom, 5/x#. The latter term represents mag- 
netization of itinerant holes. 

The applicability of the effective kinetic-exchange 
model relies implicitly on the perturbative character of 
the microscopic p — d hybridization. The level of the 
p — d hybridization over the studied doping range is illus- 
trated in Fig. 0] where we show the orbital composition of 
m T l BA- The filled symbols correspond to including spec- 
tral weights from all spd orbitals while the half-open and 
open symbols are obtained after projecting onto the d and 
sp orbitals, respectively. If no hybridization was present, 
then m % TBA projected on the d-orbitals would equal to the 
total m % TB A and the sp-orbital projected m™g A would 
vanish. In our TBA/CPA calculations, the <i-orbital pro- 
jected TOj^BA i s reduced by only 10% as compared to the 
total rniPg A and, therefore, the p — d hybridization can 
be regarded as a weak perturbation. The nearly constant 
value of the d-orbital projected rn™g A also suggests that 
the kinetic- exchange coupling parameter J p( i in the ef- 
fective spin Hamiltonian is nearly independent of doping 
over the whole range of Mn and hole densities that we 
study. 

The decrease of ttitba in Fig. [3] with increasing p/Nmu 
clearly demonstrates antifcrromagnetic p — d coupling 
over the whole range of dopings. The initial common 
slope for data corresponding to different Mn concentra- 
tions reflects the half-metallic nature of the hole system 
(only majority hole band occupied) when spin-orbit inter- 
actions are neglected. Here the hole contribution to mag- 
netization per volume is proportional to p, i.e., magneti- 
zation per Mn is proportional to p/Nmu- The change in 
the slope of ttitba at larger hole densities, which now be- 
comes Mn-density dependent, reflects population of the 
minority-spin hole band and, therefore, the additional de- 
pendence of hole magnetization on exchange splitting be- 
tween majority- and minority-hole bands. In this regime 
the hole magnetization per volume is approximately pro- 
portional to pNMn, he., a common trend for different 
Mn densities is obtained when ttitba is plotted versus p 
rather than p/Nmu- 

Within the TBA/CPA model, ttitba corresponds to 
the total magnetization per Mn measured by a SQUID 
and we therefore expect these experimental values to lie 
between 4 and 5 ixb and to increase with decreasing num- 
ber of holes in collinear (Ga,Mn)As ferromagnets. The 
XMCD data on the other hand reflect local and d-state 
projected contribution from Mn and should be compared 
with the half-open symbols in Fig. showing a negli- 
gible dependence on the hole density. These trends are 
indeed confirmed in our experiments. Before discussing 
the experimental data we refine, quantitatively, the above 
theoretical predictions for the total magnetization. For 
example, we expect that the number of minority holes at 
a given total hole density is underestimated in our TBA 



model. This is caused in part by the quantitative value 
of the exchange spin-splitting of the valence band in the 
TBA/CPA calculations which is a factor of 1.5-2 larger 
than value inferred from experiment. Also since SO cou- 
pling is not included in our TBA model, we obtain three 
majority bands that are degenerate at the T-point, in- 
stead of only two bands (heavy- and light-hole) as in the 
more realistic SO-coupled band structure. (This substan- 
tial deficiency is common to all calculations that neglect 
spin-orbit coupling.) In addition to the underestimate 
of the minority hole density, the TBA/CPA calculations 
also omit the reduction of the mean spin-density in the 
majority band caused by the SO-coupling and, therefore, 
the total TBA magnetization values are too small. In the 
following section we attempt to correct for these quan- 
titative shortcomings of the TBA calculations by taking 
the experimentally measure d 1 ^ 2 ? value of the p — d cou- 
pling constant, J p d = 54 ± 9 meV nm 3 , and by evaluat- 
ing valence band spin-splitting and SO coupling effects 
within the semiphcnomenological k • p kinetic-exchange 
modeliSiiSa The weak dependence of the TBA/CPA va- 
lence band splitting on positional disorder justifies our 
use of the virtual crystal approximation in this semiphc- 
nomenological modeling of collinear (Ga,Mn)As ground 
states. 
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FIG. 3: Mean-field total magnetization per Mn as a func- 
tion of hole density relative to the local Mn moment density. 
Results are obtained using the TBA/CPA model. 



B. Mean- field magnetization contributions from 
hole kinetic and exchange energies 

Within the semiphcnomenological virtual crystal 
model the valence band holes experience a mean-field, 
fiMF — JpdNMn(S), and the band Hamiltonian can then 
be written as, Hmf = Hkl(B) + 1imfs z , where Hkl(B) 
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FIG. 4: Integrated total and d- and pd-projected magnetiza- 
tions per Mn as a function of hole density relative to the local 
Mn moment density. See text for definition of ttl'tba- 



is the -B-dependent six-band Kohn-Luttinger Hamilto- 
nian of the GaAs host bandS 2 * 2 ^, s z is the ^-component of 
the hole spin operator, and (5*) is the mean spin polariza- 
tion of the local Mn moments. 2 W 4 At T = 0, (S) = 5/2 
and the local moment contribution to the magnetization 
per Mn is 5fiB- As emphasized above, this finding is not 
in contradiction with the smaller d-electron contribution 
to the magnetic moment in microscopic calculations. 

Because of the SO interaction, both orbital and Zee- 
man couplings of the external magnetic field have to be 
included in Hkl{B). The SO-coupling and heavy- hole - 
light-hole mixing at finite wavevectors lead to magneti- 
zations that cannot be expressed using a constant, Mn- 
and hole-density independent g-factor. Instead the ki- 
netic band energy contribution to MF magnetization per 
Mn, rrij^p, is obtained by numerically integrating over all 
occupied hole eigenstates of H MF and by finding the co- 
efficient linear in B of this kinetic energy contribution to 
the total energyiSS Results for several typical Mn dopings 
and hole densities are shown in Fig. [5] They agree quan- 
titatively with earlier calculations reported in Ref. |U 
As anticipated in Section [H] rn]^%, is negative, i.e. it 
suppresses the total magnetic moment. The magnitude 
of the hole MF magnetization per hole, \rrij^p\NMn/P: 
is smaller than 2fig due to occupation of both majority 
and minority heavy- and light-hole bands at these typi- 
cal (Ga,Mn)As hole densities (see inset of Fig. |SJ) . In this 
case, as also emphasized in Section [HI m^F is expected 
to fall into a common trend for different Mn densities 
when plotted against p. Data shown in the main panel 
of Fig. confirm this expectation and indicate a —0.2 
to QAfiB suppression of the MF moment per Mn due to 
the hole kinetic energy contribution to magnetization. 22 

To estimate the hole exchange energy contribution to 
magnetization per Mn, m^ F , we use an expression of 



the total exchange energy derived in the absence of SO- 
coupling and assuming spin-up and spin-down heavy-hole 
bands with effective mass 0.5m P , 



2 l/3 



D 4/3 



4/3, 



(6) 



where p-\(i) is the density of the majority(minority) band, 
Pl+Pl = p, and the exchange energy of the paramagnetic 
state is given by 



1/3 



p 



4/3 



(7) 



The independent hole densities can be written as 

2G T G+ 



p r (B)=pi(0)-\g jf i B Bj 
Pl (B) = Pl (0) + \ gj » B Bj\ 



2G T G^ 



(8) 



where Glw is the density of states at the Fermi energy 
of the majority (minority) band. Combining Eqs. ©-(JSJ 
and Eq. we obtain 



m ex 
171 MF 



4 2V3 2G^G* 



GT +Gi 



(p T (0) 1/3 -fi(0) 



1/3 



(9) 

In Fig. [f)| we plot rn e ^ F values calculated assuming spin- 
1/2 holes, i.e. gjj = 1. This contribution to the total 
MF magnetization is again negative, is nearly indepen- 
dent of x and p within the range of doping considered, 
and its magnitude is about a factor of 10 smaller than 
the magnitude of the term originating from the hole ki- 
netic band energy. A more realistic estimate of 
be obtained by using gjj = g/u * 3/2 in Eq. ©. Here 
g/u follows from fitting the TUjVp calculated within the 
parabolic heavy-hole band model to the full six-band nu- 
merical results in Fig. [5] The values of gfu as a function 
of x and p are plotted in Fig. [7| These are similar to 
the 5j=3/2 = 4/3 value that follows from the local atomic 
model and, therefore, m\ x IF calculated from Eq. © us- 
ing gjj = gfu * 3/2 is approximately a factor of 2 larger 
than m e ^ IF calculated assuming spin-1/2 holes. Combin- 
ing all these considerations we conclude that the zero- 
temperature magnetization per Mn in the MF kinetic- 
exchange model has a positive contribution equal to 
from the Mn local moments and a negative contribution 
from band holes which suppresses the moment per Mn 
by -5-10%. 



C. Quantum fluctuation contribution to the 
magnetization 

In Section [H] we argued that quantum fluctuation cor- 
rections to the isolated Mn d 5 +hole complex magne- 
tization should be small. Here we demonstrate that 
this conclusion also applies to the many-Mn system. 
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FIG. 5: Mean-field kinetic energy contribution to the hole 
magnetization per Mn as a function of hole density. These 
results were obtained using the six-band Kohn-Luttinger pa- 
rameterization of the valence band and the kinetic-exchange 
model. 



FIG. 7: Effective hole Lande g-factor obtained by fitting the 
kinetic energy term in hole magnetization calculated with the 
parabolic band kinetic-exchange model with j — 3/2 and 
m* = 0.5m e to the numerical results of the six-band Kohn- 
Luttinger model from Fig.|5] 
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FIG. 6: Mean-field hole-hole exchange energy contribution to 
the hole magnetization per Mn as a function of hole density. 
These results were obtained using the spin-1/2, m* = 0.5m e 
parabolic band kinetic-exchange model. 



In these calculations we use the virtual crystal kinetic- 
exchange model and assume spin-1/2 heavy-holes with 
no SO-coupling and with the parabolic band dispersion 
(to* = 0.5m e ). The many-body Hamiltonian of the 



model reads 
H = I d 3 r 



- 



ft 2 V 2 
2m* 



\J pd \J2S(Ri)-](mr-Ri) 
I 

I 



(10) 



The imaginary time path-integral formulation of 
quantum statistical physics combined with a Holstein- 
Primakoff bosonic representation for the Mn local mo- 
ments allows^ us to formally express the free energy of 
interacting local and itinerant spins in terms of a path 
integral over coherent state labels z, z: 



Z = J T>(zz) exp (-S e s[zz]) . 



(11) 



The effective action S s in Eq. (fTT|l is obtained by 
integrating out fermionic (hole) degrees of freedom in 
Eq. 110(1 . In the Gaussian fluctuation approximation 25 



5eff \zz\ 



1 \ f d 3 q 

P m */|q|<8e 



Here the inverse of the spin- wave propagator Z?(q, v m ) is 
given by 

D^ 1 (q,ii' m ) = -iv m + Ssw(B) + £ sw (q, iv m ,B) , (13) 

q c = (G^A^n) 1 ^ 3 is a Debye cutoff which ensures that 
we include the correct number of local-moment degrees 
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of freedom and 



(B) = -gsHBB 



I J, 



pd 



(PT(B)-p;(B)) (14) 



is the mean-field local moment spin-flip energy. 
The frequency-dependent self-energy E su) (q, ifl, B) in 
Eq. fT3|) is given by 



E su) (q,in,B) 



NMnJp d S 



d 3 k /( Ek - M + A(g)/2) - /( £k+q — \i — A(B)/2) 
(2tt) 3 ift + e k - e k+q + A(B) 

with e k = H 2 k 2 /2m* and 

A( J B)=iV Mn |Jprf|S'- 



(15) 



The translational and rotational invariance of our model 
implies that E su) (q, ifl) depends only of |q|. The func- 
tional integration in Eq. (|llf> can be performed exactly 
using Eq. lfT2|). 



n, |q| <q c 



i 



£ sw 

(B) +E sw (q,iv n ,B) 



(16) 



The quantum fluctuation correction to the free energy 
then reads, 



Woe 



F-F 



1 ^ Z 



(17) 



where -Zo is given by Eq. (|16[1 with the fluctuations term, 
^sw(Qi Wn, B), set to zero. 

In Fig.|Slwe plot the quantum fluctuations contribution 
to the magnetization per Mn obtained from 



asE t 



mQF 



QF 



dB 



(18) 



B=0 



As expected itiqf is small in the many-Mn system and we 
can conclude that quantum fluctuations lead to a ~1% 
suppression of the MF moment per Mn. 



IV. EXPERIMENTAL DATA 

A. Sample growth and preparation 

In the theoretical sections of this paper we eval- 
uated the zero-temperature magnetization per Mn in 
(Ga,Mn)As ferromagnets, using approximation schemes 
that would fail if the true ground state magnetiza- 
tion was highly non-collinear. We will now show that 
these theoretical results are consistent with low temper- 
ature magnetometry and XMCD experiments. A se- 
ries of (Ga,Mn)As films with Mn content varying be- 
tween 1.7-6.7% in the SQUID experiments and between 
2.2 and 8.4% in the XMCD experiments were grown by 
low-temperature molecular beam epitaxy (MBE) using 




FIG. 8: Quantum fluctuation contribution to magnetization 
per Mn as a function of hole density. These results were ob- 
tained using the spin-1/2, m* = 0.5m e parabolic band kinetic- 
exchange model. 



As2- The layer structure of the thin films consists of 25 
or 50nm (Ga,Mn)As / 50nm low temperature GaAs / 
lOOnm high temperature (580°C) GaAs / SI-GaAs(lOO) 
substrate. The growth temperature of the (Ga,Mn)As 
layer and the GaAs buffer was 180-300°C, decreasing 
with increasing Mn concentration in order to minimize 
As antisite densities while maintaining two-dimensional 
growth and preventing phase segregation. Further details 
on the growth are presented elsewhere ^SiSi 

The Mn concentrations were deduced from the in-situ 
measured Mn/Ga incident flux ratio, which was cali- 
brated using secondary ion mass spectrometry (SIMS) 
measurements on 1/im thick (Ga,Mn)As films, grown un- 
der otherwise identical conditions to the samples con- 
sidered here. A detailed comparison of the results of a 
number of different calibration techniques, presented in 
detail elsewhere^ allows us to assign an uncertainty of 
10% to the quoted total Mn doping, x. The SIMS mea- 
surements yield no information on the lattice site of the 
incorporated Mn, and we expect that Mn will be incor- 
porated either on interstitial Mni or on substitutional 
MnQa sites^S Post-growth annealing of the samples is 
performed in air at 190°C for 50-150 hours, which is 
an established procedure for removal of Mni from the 
(Ga,Mn)As layer. Curie temperatures in the as-grown 
materials are within the range of 40-80 K and in the an- 
nealed samples between 40 and 150 



B. Magnetometry 

The magnetic moment of the samples is measured in a 
SQUID magnetometer, at 5K and under a 0.3T external 
magnetic field. The external field is necessary to over- 
come in-plane anisotropy fields, so that the magnetiza- 
tion is aligned with the measurement axis of the SQUID. 
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The diamagnetic contribution from the substrate is sub- 
tracted. Measured magnetic moments normalized to the 
total Mn concentration as obtained from SIMS calibra- 
tion, msQuiD, are shown in Fig. [§] The moment de- 
creases with increasing Mn concentration, and increases 
on annealing, similar to earlier reports^ This is consis- 
tent with the anticipated formation of interstitial Mn for 
doping above ~2%, 12 given the antiferromagnetic cou- 
pling between Mni and Mnc a r^ and with breaking of 
this coupling by low-temperature annealing^^i 

In order to compare the experimental data with the 
theoretical results of previous sections we have to replot 
the measured magnetizations as a function of the density 
of uncompensated Mnc a local moments, x e ff- To do 
this we need to determine the densities of substitutional 
MnQa and interstitial Mni, x s and Xi, in our (Ga,Mn)As 
materials. Given these values, we assume^ that each 
Mni present in the system is antiferromagnetically cou- 
pled to one MnQa and that both should be excluded from 
the active Mn fraction for comparison between theory 
and experiment, i.e. that x e ff = x s — X{. 

To obtain the individual Mn impurity concentrations 
we rely on Hall effect and magnetoresistance measure- 
ments at high magnetic field (up to 16. 5T) and low tem- 
peratures (down to 0.3K), from which we evaluate the ex- 
perimental hole density p, after using a fitting procedure 
to separate normal and anomalous contributions to the 
Hall resistance^ We then assume that the single accep- 
tors MiiGa and double donors Mni are the only impuri- 
ties that contribute to p, i.e., p = (4/af c )(x s — 2xj). From 
this expression and from the total Mn concentration ob- 
tained by SIMS calibration (x — x s +Xi) we can estimate 
x e ff for both as-grown and annealed samples. A detailed 
discussion of the uncertainties associated with this proce- 
dure is given elsewhere. 12 The magnetic moment per ef- 
fective Mn moment density, m'sQUiDi 1S snown i n Fig. llQI 
as a function of p/N^ where N%££ — (4x e ff/af c ). In 
agreement with the predictions of the theory section, 
m SQUW ^ aus w ithin the range 4-5/ie for all samples 
studied. Furthermore, although there is appreciable scat- 
ter, it can be seen that samples with lower hole densities 
tend to show higher rn e gQ UID , consistent with a negative 
contribution to magnetization from antiferromagnetically 
coupled band holes. 




FIG. 9: SQUID magnetization per nominal total Mn den- 
sity in as-grown (open symbols) and annealed (filled symbols) 
(Ga,Mn)As materials plotted as a function of the nominal Mn 
doping. 
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C. X-ray Magnetic Circular Dichroism 

XMCD measurements were performed using (99 ± 1)% 
circular polarized x-rays from beamline ID8 at the Eu- 
ropean Synchrotron Radiation Facility. The samples are 
briefly etched in concentrated HC1 prior to the measure- 
ments in order to remove Mn oxide rich surface layers, 
which may obscure the signal from the (Ga,Mn)As due to 
the relatively short probing depth of the measurement^ 
After etching, total electron yield and fluorescent yield 
measurements are in quantitative agreement, indicating 



FIG. 10: SQUID magnetization per effective density of un- 
compensated MnGa local moments in as-grown (open sym- 
bols) and annealed (filled symbols) (Ga,Mn) As materials plot- 
ted as a function of hole density per effective density of un- 
compensated MnGa local moments. 



a uniform distribution of Mn. 

Fig. II II shows Mn 1^2 x-ray absorption spectra for an 
annealed (Ga,Mn)As sample with x = 8.4%, for paral- 
lel and antiparallcl orientations of the external magnetic 
field and the x-ray helicity vector. The sample temper- 
ature is 6K, and the external magnetic field is ±1T, ap- 
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plied perpendicular to the surface of the sample. A very 
large change in the absorption is observed on reversing 
the external field, with an asymmetry (difference to sum 
ratio) of up to 55% at the L3 peak. absorption cor- 
responds to transitions from the 2p core states to the un- 
filled 3d states, so the Mn 1,3,2 spectra gives direct infor- 
mation on the polarization of the Mn 3d band. Applying 
the XMCD sum rules to the spectra allows quantitative 
and separate determination of the Mn 3c? ground state 
orbital and spin magnetic moments^ The moments are 
obtained on a per atom basis, without requiring separate 
measurement of the Mn concentration, by normalizing to 
summed absorption signal. There are, however, inherent 
uncertainties in the application of the sum rules, in par- 
ticular due to mixing of the 2p 3 / 2 an d 2p%/2 states which 
prevents the separate integration over each of the spin- 
orbit split core levels. Comparison of calculated spectra 
with their corresponding ground state moments reveals 
that a correction factor of 1.47 is required for the spin 
moment to account for this mixing. The Mn 3d mo- 
ments obtained from XMCD are shown in Table I, for 
two annealed samples with low and high Mn doping. In 
both cases, magnetic moments of around 4.5/xs are ob- 
tained, in agreement with the SQUID results. Moreover, 
the measured moment and the property that it is inde- 
pendent of Mn doping is in very good agreement with the 
calculated d-projected magnetic moment shown in the in- 
set of Fig. 0] We note here that the calculations shown 
in Fig. 01 account only for the spin angular momentum 
contribution to magnetization since SO-coupling effects 
were neglected in the TBA/CPA calculations. 
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0.15 
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0.16 



4.5 



TABLE I: Mn 3d moments obtained from XMCD and de- 
composed into the spin and orbital contributions in annealed 
samples with nominal Mn doping 2.2 and 8.4%. 



V. CONCLUSIONS 

We report on a combined theoretical and experimental 
analysis of the spontaneous magnetization in (Ga,Mn)As 
diluted magnetic semiconductor ferromagnets. We find 
that the thermodynamic magnetization is dominated by 
a large local moment contribution of 5/j.b from nearly 
collincar substitutional Mn atoms. Evaluation of the 
smaller magnetization contribution from the valence 
band system that couples the local moments together in- 
volves a number of subtleties. In this paper we included 
hole-hole exchange interactions and we also accounted for 
spin-orbit coupling which means that no valence band 
orbital is completely spin-polarized and which substan- 
tially changes the overall electronic structure. Quantum 
fluctuations of the band and local moment orientations 
also play a role because of the antiferromagnetic inter- 
action between band and local moment spins. The end 
result of all these corrections is a magnetization per mag- 
netically active Mn ion that is suppressed from 5/ib by 
~ 5 — 10%. Comparison with experimental data can be 
made reliably only after accounting for the formation 
of interstitial Mn complexes during the MBE growth, 
and for their subsequent removal by post-growth anneal- 
ing. Once these corrections have been applied, we find, 
within the experimental error bars, agreement between 
theory and experiment. The interpretation of XMCD 
magnetization measurements, which capture only the d- 
electron contribution, requires a recognition of the hy- 
bridized p-d character of both local moment and band- 
electron contributions to the magnetization. Comparison 
of these measurements with TBA/CPA calculations pro- 
vides experimental support for the applicability of the 
kinetic-exchange model in (Ga,Mn)As ferromagnets. Fi- 
nally, our combined theoretical and experimental work 
demonstrates that non-collinearity does not play a sig- 
nificant role in the magnetization of high-quality metallic 
(Ga,Mn)As ferromagnets. 
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